Wobbling motion in the multi-bands crossing region 
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Abstract 

The backbending in the A ~ 180 mass region is expected to be caused by multi- 
bands crossing between low- if (g- and s-bands) and high- if bands (K n = 8 + or 
10 + ). We analyze a mechanism of coupling of these bands in terms of a dynamical 
treatment for nuclear rotations, i.e., the wobbling motion. The wobbling states 
are produced through the generator coordinate method after angular momentum 
projection, in which the intrinsic states are constructed through the 2d-cranked 
HFB calculations. 
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A new interest in the "backbending" phenomenon has arisen through a series of 
experiments in 1990s' [1-4] for nuclei in the A ~ 180 region. It is expected that 
the backbending in this mass region is caused by the three-bands crossing, i.e., 
the ground (g-), super (s-) and "tilt" (t-) bands [4,5], unlike the backbending 
caused by the two-bands crossing of the g- and s-bands [6] in light rare-earth 
nuclei. The t-band was recently proposed as "tilted rotational band" [7]. In 
nuclei in the A ~ 180 region, the Fermi level lies in states of high-f2 (Q is a 
projection of single-particle angular momentum onto the symmetry axis of a 
nucleus). As a result, low-f2 states associated with the rot at ion- alignment are 
mostly occupied. Therefore, it is natural to expect a new type of excited bands 
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involving the high-f2 states. Hereafter we refer to the bands as high-.fr bands. 
Such bands as K n = 8 + (10 + ) are actually observed in the vicinity of the yrast 
line □ in 180 W and 182 0s [1,2] ( 182 W and 184 0s [3,4]). It is characteristic that 
these bands have inter-band E2 transitions towards the yrast band having 
low- A" configurations, namely these transitions violate the A"-selection rule to 
a great extent (A A" ~ K). Walker et al. considered that these transitions are 
related to a new type of the backbending mechanism involving tilting degrees 
of freedom in the rotation axis [1] . The nature of the high- A' bands changes at 
lower and higher spin than that in the backbending region (J c ~ lAh). At I < 
I c , the bands form a regular AI=1 sequence, while at / > I c the bands split 
into two AI=2 sequences. After the splitting, the odd-spin sequence becomes 
energetically lower than the even-spin one, although the high-A" bands begin 
with even spin (I = 8h or lOti). This type of splitting is conventionally called 
"signature inversion". These features in the high- A' bands suggest a certain 
strength of inter-band interaction between the low- and high- A' bands. 

There are, at least, two modes that couple states having intrinsic structures 
different in the A'-quantum number. One is the wobbling motion that we 
present in this letter and the other is the triaxial deformation (i.e., sizeable 7 
deformation). In our recent calculations based on the two-dimensional cranked 
HFB (2d-CHFB) [8], the intrinsic states of 182 . 184 0s have (almost stable) pro- 
late shape (7 < 3°; (3 ~ 0.275) in the backbending region (I ~ I c ). Also, the 
calculations show that tilting and 7 degrees of freedom are almost decoupled in 
the region, while the coupling is seen in the very high-spin region like / ~ 36h. 
In this letter, we focus on the wobbling motion decoupled from the 7-degree 
of freedom in I ~ I c . 

The conventional cranking model, which is restricted to the rotation about the 
principal axis of the mass quadrupole deformation (PAR), can be extended to 
the rotation about an axis tilted away from the principal axis (TAR). PAR 
simulates low- A" bands such as the g- or s-bands, while TAR describes high- A" 
bands . In the present paper, both the PAR and TAR are calculated self- 
consistently through the 2d-CHFB method. We regard the wobbling motion 
as a dynamical mode to couple the low- and high- A' states in a quantal and 
microscopic approach. We treat the motion in terms of the generator coordi- 
nate method after the angular momentum projection (GCM after AMP), in 
which tilt angle of the rotating axis is chosen to be the generator coordinate. 

In our previous paper [11], we presented calculations based on the wobbling 
model (GCM based on 2d-CHFB wave functions) without angular momentum 
projection and encountered a difficulty in solving the Hill- Wheeler equation: 
a problem concerned with convergence in the eigenvalues with respect to the 



"yrast line" means a sequence of the lowest states in energy for given angular 
momentum. 
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"cut-off" dimension, which we will discuss below. We tried to improve the 
method through the constrained Hill- Wheeler equation, but could not obtain 
satisfactory results [12]. All the difficulties seem to come from the broken sym- 
metries by the mean field approximation. (Angular momentum and numbers 
of nucleons are not conserved in deformed HFB states.) In many cases, the 
number constraints may be good enough to produce the reasonable cranked 
HFB states of stable nuclei [13,14]. We consider that, in a study of high-spin 
states, restoration of the rotational symmetry is more important than the 
gauge invariance, as the first attempt to circumvent the difficulties. 

The 2d-CHFB states are wave packets of eigenstates in angular momentum 
(and number) and have broad width around the constrained value J [10,14]. 
Thus, the GCM states also include substantial components of undesirable 
angular momenta far from I — J. In order to eliminate them, we perform 
exact angular momentum projection (AMP) on 2d-CHFB states. Norm and 
energy overlap kernels are obtained through the formula [15]. But the well- 
known ambiguity of the phase by n inheres in the formula for the norm overlap 
kernels. We employ the analytic continuation method to determine the proper 
branches. [16]. 

The parameter set of the pairing+Q-Q model Hamiltonian employed here and 
details of the self-consistent 2d-CHFB method are found in Ref.[8]. In the 
3d-cranking model, constraints on three components of angular momentum 
Jit = (<E> | Jfc | $) are specified by the tilting angle (9, <fi) and J through Ji = 
J cos 9 cos <fi, J 2 = J cos 9 sin and J 3 = J sin 9. Here, a cranked HFB state is 
solved for 182 0s and denoted as J)). The latitude (9) and longitude (0) 

angles are introduced so as to measure the deviation of the rotating axis from 
the x-axis (PAR). 

Now, we set up a wobbling state, or GCM state (I 1 ?)), which is described as 
the superposition of projected 3d-cranked HFB states with different tilting 
angle: 

l do 00 

I*m>= E fd9jd^f I K mPL K \Ho<P;J)), (i) 

K= - I -e -<p 



where Plj K is the angular momentum projection operator [14]. The function 
fki^i 0) i s called the generator wave function. The azimuthal angle O is taken 
to be zero in the present calculations and is omitted hereafter. As for 9$ on the 
integration limit, we will explain later. J is fixed to be 14h, and we also see 
later some properties of the CHFB states with a constraint J = lAh briefly. 
The GCM state is obtained by the variational principle with respect to fx(9), 
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leading to the so-called Hill- Wheeler equation: 

E / de' (H^e, e>) - EfN^e, e>)) = o. (2) 



The norm and energy overlap matrix-kernel N I KK ,{9,9') and H I KK ,{9,9') are 
respectively defined as follows, 



(3) 



We solve Eq.(2) in two steps [11]: (1) Diagonalize the norm overlap matrix 
kernel like, 

E / dffN^e^x 1 *'^) = ^(O), (4) 



to obtain eigenvalues and eigenfunctions Xn K {@), m which n is an index 
specifying v J n in order of magnitude. The eigenvalues are never negative owing 
to properties of the norm overlap matrix-kernels; (2) Define energy matrices 
on the base of orthonormal set | nIM) = J2k I% Mx™ C ( d ) '/ ' \Rn p MK \^(°)) 



(5) 



Kn< = E // d9d9'^^H I K A0,0')^X 

KK' JJ \IK sJv l n , 



Eventually, the Hill- Wheeler equation gets transformed into an eigenvalue 
equation: J2 n > 7~Lnn' g n > = E 1 g n . Integrations in Eq.(5) are achieved numeri- 
cally by an approximation of discretization. In the present study, five angu- 
lar points are chosen (| <&(«/, 9));9 = 0°,±7° and ±20°) under a condition: 
< 0.9, up to |^ m ax| = 20°. There are two reasons why we take 
only angles up to 20°. The first reason comes from a physical argument: ac- 
cording to the energy curve for the CHFB states with J = lATi (See the curve 
for J = 15h in Fig. 4(a) in Ref.[8]), the first valley around 9 = 0° has an edge 
at 9 — 20°. Realization of the wobbling motion around PAR is thus achieved 
approximately by the superposition of the CHFB states corresponding to this 
first valley. Although we should attempt to examine a global wobbling mo- 
tion from 9 = —90° to 9 — +90°, the present framework is satisfactory to 
study the mechanism of mixing between high- and low-i^ states. The second 
reason comes from numerical difficulty to calculate the norm kernels at larger 
tilt angles which include quite high frequency modes. For a precise treatment 
of overlap kernels we need much finer meshes in the Euler angles, but this 
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requirement in numerical calculations cannot be feasible within an available 
computation time. 

The net dimension of the Hill- Wheeler equation turns out to be n max = 
{5 x (21 + 1)} for a given angular momentum (J). Fig.l shows the spectrum 
of each band (the yrast, yrare, and second yrare band ) as a function of n cut . 
A "plateau" indicates that the solutions (i.e., energy eigenvalues) are almost 
independent of the cut-off dimension, n cut . This result implies a great advan- 
tage of AMP in comparison with our previous GCM calculations. From the 
results, we truncate at n CVLt = 11 for all states except the I — 17 yrare state, 
in which we truncate at n — 10. The odd-/ states in the first yrare band form 
"slopes", except narrow plateaus formed around n cnt = 10. The reason for the 
appearance of the "slopes" may come from mixing of the states with the differ- 
ent particle number, i.e., "contaminations in number", for we have not carried 
out number projection in this work. Though the "contaminations" make our 
discussion dubious, the dual projection in number and angular momentum is 
nothing but a heavy calculation consuming a great deal of computer time. 
At present, we leave an examination on the effect of number projection as a 
future study. 

Fig. 2 shows the resultant energy spectrum in the band crossing region. The 
signature splitting and inversion take place in the / > 15h region. The phase of 
the staggering is properly reproduced comparing with the experiment: odd-/ 
states are energetically lower than even-/ states. However, the amount of the 
splitting is much larger than experimental values. Our result is roughly 500 
keV while the experiment is about 100 keV. There are several reasons for this 
discrepancy. There may be some argument for a choice of the pairing-plus- 
quadrupole force, but a more important factor is a set of parameters for the 
effective force. The parameters should be adjusted for calculations of the GCM 
after AMP states, though this procedure is again difficult because of the limit 
of the present computation resources. Thus, we have to keep using the same 
parameters as in the HFB calculations. Other possibility for the discrepancy 
is that we take a fixed J value (J = lAh) in the CHFB states. For example, 
in this study, the projected states with / = 17h are projected out from the 
CHFB state with J = 14h. To make a better projected states, it can be better 
to project out from the CHFB states with J = 

It is satisfactory that the A/=l band sandwiched between the two A/=2 bands 
is reproduced^ in Fig.2. We are tempted to say this band structure is pro- 
duced via multi-bands crossing among g, s, and high-K bands. However, we 
should note that the higher A/=2 band does not have a g-band character but 
a s-band one. This result can be explained as below: we have fixed the an- 



3 "yrare band" means a sequence of the first excited states from the yrast states. 

4 We pick up the lowest three levels out of 11 levels in order to draw the Fig.2. 
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gular momentum constraint on J = 14ft, and the gap energy for neutron in 
this range of angular momentum is substantially reduced due to the rotation- 
alignment (Refer to the result in Ref.[8]). Therefore, our generating states 
have a good amount of the rotation-aligned component, i.e. a strong character 
of the s-band. Consequently, the upper Ai=2 band is of the s-band character 
instead of the g-band. 

Backbending is not well reproduced in the present calculations. We can sup- 
pose that this deficiency is also due to the lack of the g-band components 
in our calculations because the g-band can contribute to up- or back- bend- 
ing through its crossing with the s-band. To make a complete analysis of the 
three-bands crossing in 182 0s, we need to use, at least, two different CHFB 
states: one with low J and the other with high J, and perform the GCM after 
AMP calculations around at each angular momentum. However, for a practi- 
cal reason, the present paper is dedicated to an intensive study of the mixing 
of low- and high-i^ states through the wobbling motion, which corresponds to 
the band crossing of the s-band with the t-band. This situation can happen 
just after the g-s crossing. As for the g-t crossing, we will discuss later. 

Let us learn the obtained K distributions of the wave functions in the context 
of the interaction between the lowest two bands. For this purpose, we calculate 
overlaps between the wobbling (GCM) state and each spin-projected intrinsic 
state ( i.e., 2d-CHFB states): (<E»(J, 0)\Pmk\^m) = 9k(°) [ 14 L evaluated by 
using the formula, 

"cut 

^) = EVJw»rt- (6) 

n=l 

Fig. 3 shows the quantity |(?i-(#)| 2 with respect to K for each 9 value. In our 
calculations, the / = 14ft and / = 16ft states correspond to the states before 
and after the crossing, respectively. In / = 14ft, the yrast state shows a typical 
feature of the s-band: the major component is K = and some fluctuations 
around it, i.e., ~ lft or 2ft. This feature implies that a good amount of the 
rotation-aligned components are induced in the 2d-CHFB states at J = 14ft, . 
In turn, the yrare state has two large peaks at K — ±8ft, and one small peak 
at K = Oft. These features indicate a typical high-fC band. Most of the high- 
K components are brought about by TAR states (\9\ = 20°). The existence 
of the K — component implies the inter-band interaction with the yrast 
band, though the interaction is weak due to the large difference in energy 
[8]. In the crossing region, one can expect that the mixing is supposed to be 
stronger as the two bands come closer. In / = 16ft, we can find that both of 
the yrast and yrare states have three peaks at ±8ft and K — 0, indicating 
that the bands are mixed with each other after the crossing. In particular, the 
yrare states have the three comparable peaks at low- and high-fT components. 
In these perturbed states, the high-.K" components come mainly from TAR 
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states while the low-ZT components are brought by PAR states. Hence, these 
states are interpreted as the wobbling states, namely, dynamically rotating 
states coupling the low-ZT PAR and high-ZT TAR states. On the other hand, 
odd-/ states in the high-Zv" band are much less perturbed according to the 
calculations. This is because the odd-/ members are energetically isolated and 
have no partners lying nearby to couple ( PAR has only even-/ states). 

This fact draws an interpretation for the signature inversion in the t-bands: 
The inter-band interaction between the s- and t-bands pushes up (down) the 
even-/ states in the t-band (s-band). (Note that the s-band is the yrast after 
the crossing with the g-band.) This downward shift of the yrast band may 
enhance backbending, though it is not clearly seen in the present calculation. 
The interaction between the g- and t-bands is quite weak because the g-band 
is approximately the pure K — state with small admixture of high-ZT com- 
ponents. As a consequence, the backbending in the A ~ 180 region occurs not 
only due to the rotation-alignment mechanism dominating in the rare-earth, 
but also due to the mechanism of the wobbling motion, i.e., the low- and high- 
K interaction. On the other hand, the odd-/ states are unperturbed and serve 
as reference states for the interaction energy. 

In summary, we have investigated a mechanism of backbending and the signa- 
ture inversion in 182 0s, by means of GCM after AMP, on the 2d-cranked HFB 
states. With this method, we have qualitatively reproduced the main feature 
of the level structure showing the signature inversion in the high-ZT band ly- 
ing very closely to the yrast band. We interpret this result from a point of 
view of an inter-band interaction between the low-ZT and high-ZT bands. We 
have shown that the perturbed states have characters of the wobbling mo- 
tion, that is, dynamical mode coupling low-ZT PAR and high-ZT TAR states. 
In terms of this wobbling model, we have discussed an enhancement of the 
backbending in the A ~ 180 region, in which the typical rotation-alignment is 
somewhat suppressed due to the location of the Fermi level. It is also shown 
that AMP brings about a great advantage in stabilization of eigenvalues in the 
Hill- Wheeler equation. At the same time, we realize that number projection 
turns out be important since it may eliminate a "slope" owing to the "number 
contamination" , as we expand the GCM states with the orthonormal basis. 

We would like to thank Dr. N. Tajima for discussions and suggestions on this 
work. Most of the computations have been done using the Vector Parallel Pro- 
cessor, Fujitsu VPP500/28 at RIKEN which are also gratefully acknowledged. 
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Figure Captions 



Fig.l : A relation between the cut-off dimensions (n cut ) and the energy eigen- 
values (E / (n cut )) , in the Hill- Wheeler equation. The "plateaus" of the graphs 
imply the possibility of "safety cut-off" within the plateaus. In the present 
study we choose n cut = 10 for the first yrare I = 17h state and n cut = 11 for 
the other states. 

Fig. 2 : Energy spectrum calculated through the GCM after AMP. The angu- 
lar momentum constraint for the intrinsic state is fixed to be J = 14h. From 
the analysis of the generator wave functions, the yrast and the second yrare 
bands (AI=2 -bands) are of the low-.fr characters, while the first yrare band 
(Ai=l -band) is of the high-if character. The signature inversion is seen at 
the region I > 15h. 

Fig. 3 : Graphs for the overlaps between the GCM state and the spin-projected 
intrinsic state: g^iO) = (^(J, 9)\Pmk\^m)- (See texts.) The left panels are for 
the GCM states before the band crossing (I = I4h), and the right panels are 
for the states after the crossing (/ = 16h). The lower and upper panels show 
the yrast and first yrare states, respectively. In each panel, the graphs are 
separately drawn for each tilt angle: 8 = 0°, ±7°, and ±20°. 
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Stability of energy eigenvalues in the Hill- Wheeler equation 
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Yrast of the GCM states; 1=14 
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Yrare of the GCM states; 1=16 
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